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■^sf ■ Abstract 

O 
O 

Cosmic microwave background (CMB) experiments generally infer a temper- 
ed ', ature fluctuation from a measured intensity fluctuation through the first term 
in the Taylor expansion of the Planck function, the relation between the inten- 
ts. ' sity in a given frequency and the temperature. However, with the forthcom- 
ing Planck satellite, and perhaps even with the Microwave Anisotropy Probe, 
the CMB-dipole amplitude will be large enough to warrant inclusion of the 
, next higher order term. To quadratic order in the dipole amplitude, there is 
an intensity quadrupole induced by the dipole with a frequency dependence 
given by the second derivative of the Planck function. The Planck satellite 
^ should be able to detect this dipole-induced intensity quadrupole and distin- 
guish it through its frequency depdendence from the intrinsic CMB tempera- 
ture and foreground quadrupoles. This higher-order effect provides a robust 
pre-determined target that may provide tests of Planck's and MAP's large- 
angle-fluctuation measurements and of their techniques for multi-frequency 

?-h ' 

» ■ foreground subtraction. 
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I. INTRODUCTION 



The primary aim of NASA's recently launched Microwave Anisotropy Probe (MAP) |IJ 
and the European Space Agency's Planck satellite will be to measure small-scale fluctua- 
tions in the cosmic microwave background (CMB) in order to better determine cosmological 
parameters P,|J and test inflationary cosmology ||, and thus improve upon the already re- 
markable results from recent balloon and ground experiments || . However, the most salient 
feature in the CMB is not this small-scale structure; it is the much more prominent dipole 
of amplitude AT = 3.365 ± 0.027 mK. 

The simplest explanation for the dipole is a Local Group velocity v = 627 ± 22 km sec -1 
toward (I, b) = (276°±3°, 33°±3°) with respect to the CMB rest frame. Some measurements 
of galaxy velocities on large scales suggest independently, although far less precisely, a similar 
value 0. It has thus been generally accepted that the CMB dipole is due to our peculiar 
velocity, but theorists have occasionally speculated that the dipole might be due, at least in 
part, to an intrinsic temperature fluctuation [[|. 

If we had sufficiently precise multi-frequency observations to measure the frequency spec- 
trum at each point on the sky, it would alway be consistent with a blackbody spectrum, 
although with an angle- dependent temperature. In practice, the observation frequencies and 
sensitivities are limited. The temperature at each point on the sky is thus determined by 
measuring the intensity fluctuation, and then converting it to a temperature fluctuation, 
generally assuming that the temperature fluctuation is sufficiently small that the intensity 
fluctuation can be related to the temperature fluctuation by the first term in a Taylor ex- 
pansion of the Planck function, the relation between the temperature and intensity. With 
prior experiments, the temperature fluctuation has always been sufficiently small that this 
approximation has been warranted. 

In this paper we point out that with the observed dipole amplitude and Planck's (and 
possibly MAP's) sensitivity, this will no longer be a good approximation — the second-order 
term in the Taylor expansion will be non-neglible. To quadratic order in the dipole am- 
plitude, there is an intensity quadrupole with a frequency dependence given by the second 
derivative of a Planck function. Since this frequency dependence differs from that of the 
first-order term, this dipole-induced intensity quadrupole can be disentangled from the tem- 
perature quadrupole, even after taking into account a number of foregrounds, as we show 
below. This higher-order effect will provide a robust pre-determined signal that Planck and 
possibly MAP should be able to detect. It can be used to calibrate large-angle- fluctuation 
measurements and. or benchmark multi-frequency foreground-subtraction techniques. It may 
thus be a useful addition to the astrophysical point sources and annual modulation of the 
dipole that have until now served as calibration sources for CMB experiments. 

In the next Section we discuss the dipole-induced intensity quadrupole, which we refer 
to simply as the "dipole quadrupole" and abbreviate DQ. In Section III, we then show that 
the DQ should be detectable by Planck and distinguishable from the CMB quadrupole, 
even after subtracting several foregrounds. Our analysis in Sections II and III assumes that 
the dipole is due entirely to the velocity of the solar system with respect to the CMB rest 
frame. To indicate the magnitude and detectability of of the DQ, we evaluate how well 
the magnitude and orientation of the peculiar velocity could be determined — assuming the 
dipole were due to a peculiar velocity — from the DQ alone. In Section IV we show that the 
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same DQ arises even if the CMB dipole is an intrinsic temperature fluctuation, rather than 
the result of a peculiar velocity. Section V provides some closing remarks. 



II. THE DIPOLE-INDUCED INTENSITY QUADRUPOLE 



The specific intensity of a blackbody of temperature T CMB is 



C = c- 



x 



>3 



(1) 



where C = 2(A; B T CMB ) 3 (/ic)~ 2 is a constant, x' = /i^'(A; b T CM b) _1 , v' is the photon frequency, 
and h, kg, and c are, respectively, the Planck and Boltzmann constants and the speed of 
light. The photon frequency measured by an observer moving with velocity v relative to the 
blackbody is v = 71^(1 + where (3 = v/c, 7 = (1 — /5 2 ) -1 / 2 , and // is the cosine of the 
angle between the velocity and the photon direction. The specific intensity transforms as 
I' u i/v' z = Iv/v 3 , and thus the observer- frame specific intensity is 



C 



x 



g3!7(l+/3^) 



(2) 



where x = /w(&bTcmb) 1 - Expanding in /3, 



h = C- 



x 



- { [l + 0(f3 2 )] - f(x)/3n + f(x)g(x)(3 2 (ii 2 - 1/3) + •••}. 



(3) 



where f(x) = xe x (e x — 1) and g(x) = (x/2) coth(a;/2). 

The term linear in (3 is the dipole, with the appropriate frequency spectrum f(x) for the 
dipole [10], which is also more generally the frequency spectrum for a thermal fluctuation. 
However, to order f3 2 , there is another term, the DQ, with a frequency spectrum that differs 
from that from small thermal fluctuations. With COBE-DMR's sensitivity and frequencies 
(in the Rayleigh- Jeans regime), the frequency dependence was negligible. However, with 
the improved sensitivity and expanded frequency range of new satellite experiments, the 
frequency dependence will become detectable. This frequency dependence was noted in 
Refs. , but seems to have since been overlooked. Refs. []T3[ identified the effect (sometimes 
referred to as the "kinematic quadrupole" ) in calculations of CMB polarization induced by 
reionization, but they did not consider the CMB dipole. Ref. [nj] discusses peculiar-velocity 
effects on the CMB, but only on small-scale fluctuations. 



III. DETECTABILITY OF THE DQ 

To assess the detectability of the DQ, we determine how well the amplitude and direction 
of the solar-system velocity can be determined with MAP and Planck under the assumption 
that the velocity is to be reconstructed entirely from the DQ. 

Since the frequency dependences of the CMB quadrupole and the DQ differ, they can be 
distinguished if the quadrupole moments of the intensity are determined in several frequen- 
cies. However, the CMB quadrupole and the DQ will also have to be distinguished from 
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Frequency v (GHz) 

FIG. 1. The frequency dependence of the cosmological (solid; black), kinematic (short-dash; 
blue), dust (long-dash; magenta), synchrotron (dot-short-dash; cyan), bremsstrahlung 
(dot-long-dash; red), FIRB (short-dash-long-dash; green), and Sunyaev-Zeldovich quadrupole mo- 
ments, all referenced to the frequency dependence of the cosmological quadrupole. For the Sun- 
yaev-Zeldovich effect, we plot the absolute value as the frequency dependence changes sign near 
v = 220 GHz. The amplitudes are arbitrary, as our analysis assumes that they will be determined 
by the data. 
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the quadrupole moments due to unsubtracted, or imperfectly substracted, foregrounds, each 
with its own frequency dependence. 

In what follows we discuss the brightness temperature (which is proportional to the inten- 
sity) rather than the intensity so that we work with units that are familiar in CMB studies. 
If the brightness temperature T v (n) is measured in several frequencies v as a function of 
position n over the entire sky, then the five (for m = —2,-1,0,1,2) quadrupole compo- 
nents T m>u can be constructed from T m ^ = J Y2 m (h)T u (h)dh, where Y"; m (n) are spherical 
harmonics. The possible contributions to these T m ^ from several unsubtracted foregrounds 
(dust, synchrotron radiation, bremsstrahlung, the far infrared background, and the Sunyaev- 
Zeldovich effect), which will have frequency dependences shown in Fig. [I], must also be con- 
sidered. The quadrupole moment at each frequency is then the sum of the contributions from 
all of these sources: T m>u = T^/j(i/), where % = {CMB, DQ, dust, synch, brem, firb, SZ}, 
and fiiy) is the frequency dependence of source i, referenced to the CMB frequency depen- 
dence. Thus, /cmbM = 1; /dqM = (z/2) coth(x/2); / duBt (i/) = x 2 e"(e* - l) 2 (e*« - l)" 2 
(where Xj = /i^(fcsT dust ) _1 and we have taken a thermal dust spectrum with T dust = 20 K 
and emissivity index a = 2); f sync h_(v) = x~ 4 ' s (e x — l) 2 e~ x ; /brem(^) = x~ 4Ae (e x — l) 2 e~ x ; 
ifirb(^) = x 1 ' 64 e Xf ~ x (e x — l) 2 (e Xf — 1)~ 2 (where Xf = hviksT^)^ 1 and we have taken a 
thermal FIRB spectrum with = 18.5 K and emissivity index a = 0.64 as a fit to data 
15|); and /sz( I/ ) = 2/cmb( z/ ) — /dq(^)- For a review of foreground properties see Ref. |16| . 



For each of the five components m, we will have five (for the MAP frequencies of 22, 30, 
40, 60, and 90 GHz) or nine (for the Planck frequencies of 30, 44, 65, 100, 143, 217, 353, 
545, and 857 GHz) data points T TO) „ that we model with a vector of parameters s» = 
again for i = {CMB, DQ, dust, synch, brem, firb, SZ}. Following Ref. ||, the standard error 
(la) to the quadrupole component will be cr^ Q = [ar ly ' 2 ]DQ,DQ> where 



BT 8T 

dsi dsj 



E «?MM (4) 



(7 — U 

" m,u v " m,v 



is the covariance (Fisher) matrix. Here, a 2 ^ v is the variance to the quadrupole from instru- 
mental noise at frequency u; we use 4-year values for MAP from Ref. [|IJ and 2-year values 
for Planck from Ref. 0. 

Once the DQ components have been determined, the peculiar-velocity compo- 

nents can be found. Expanding the angular dependence, (/i 2 — 1/3), where \i is the co- 
sine of the angle between the velocity and the photon direction, in spherical harmonics 
gives us (T^ /T CM b) = -^An/45((3 2 + (3 2 - 2(3 2 ), v^(R^i/?cmb) = -Jlfa/lEpA, 
^(ImT^/TcMB) = - yW/15 /y^, V2(ReT%%/T CMB ) = -^4tt/15(/3 2 - /3 2 ), and 
\/2(ImT°? 2 /TcMB) = —ijl6ir/15f3 x /3 y . Let us suppose that these moments have been 

measured precisely. Then the equations for T°= (= a), ReT°= 2 (= b), and ImT^S 2 (= c) 
can be inverted (after introducing the shorthand a, b, and c) to give us the compo- 
nents (3 2 = v /l5(16vr)- 1 ( v / & 2 + c 2 + b), [3 2 = JlhilQu^^b 2 + c 2 - b), and (3 2 = 



15(167r) _1 ("\/fr 2 + c 2 — Voa). Since we have now determined the squares of three com- 
ponents of the velocity, we are still left with a residual eightfold degeneracy in the velocity. 
However, the signs of the components ReT^^ and ImT^^, which are proportional, respec- 
tively, to (3 X (3 Z and /3 y p z , can be used to break two of these degeneracies, leaving us with only 
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a twofold degeneracy, the sign of the velocity, which is undetermined by the quadrupole. 

When the measurements are done, we will then want to check whether the peculiar 
velocity induced from the DQ is consistent with that inferred from the dipole itself. If 
the data analysis and our understanding of the instrument are reliable, then the two ve- 
locities should be consistent. To quantify the degree to which MAP and Planck can test 
this consistency we calculate, assuming v = 0, the expectation value of x 2 f° r the hy- 
pothesis that v = 627 km sec -1 . If, e.g., x 2 = 9 then the v — hypothesis can be 
distinguished from the v = 627 km sec -1 hypothesis at 3cr. We find, for arbitrary v, 
X 2 = 9(4.2 fiK/a^) 2 t>/(627 km sec -1 ) . From this x 2 we a l so derive the minimum v 
distinguishable at 3cx from v = to be v m - m = 305 km sec" 1 (a^/ziK) 1 / 2 . 

As another indicator of the precision with which the DQ can be reconstructed, we now 
determine the precision with which the magnitude and orientation of the velocity can be 
inferred from the DQ alone. There will be a measurement error of to each of the T®®, 
as described above, and the three velocity components, v x , v y , and v z , will be fit to all five 
of these quadrupole components. To estimate the errors on these components we once again 
calculate a Fisher matrix, this time for these five data points and the parameters v x , v y , 
and v z . We then invert the Fisher matrix and choose the z axis along the inferred velocity 
so that v = v z . We then find that the measurement error (at 3<r) to the velocity is Av = 



3y 45(647r) -1 a^^ 2 (t> 2 T CM B) \ where v is the best-fit velocity. Under the null hypothesis 
that the peculiar velocity is that inferred from the dipole, we find Av = 225(er^//xK) 
km sec -1 . The Fisher-matrix analysis tells us that the measurement errors to the x and y 



components (those perpendicular to the best-fit velocity) are each y4/3 times that for the z 
component. Thus, under the same null hypothesis, the (3<r) error to the velocity orientation 
will be 24°(a£%K). 

Table | shows results of our calculations of x 2 and Av assuming a variety of combinations 
of components will be fit to the data. If the dust contribution to MAP must be determined 
from the MAP data itself, then the smallest detectable peculiar velocity is ~ 2100 km sec -1 , 
too big to be interesting. This is because the highest-frequency channel becomes a dust 
monitor, and the remaining MAP channels at ^ 60 GHz do not provide enough leverage 
to disentangle the DQ and the CMB quadrupole. However, if the dust quadrupole can 
be determined precisely from other observations, then MAP might be able to isolate the 
DQ to better than 7a. The synchrotron and bremsstrahlung foregrounds are strongest 
at low frequencies, where the DQ and CMB frequency dependences are similar. Thus, if 
synchrotron and bremsstrahlung emission are included in the analysis, the lowest-frequency 
channels become foreground monitors and the ability to separate the DQ is not degraded 
significantly. 

Because of the improved detector sensitivity, and especially the broadened frequency 
coverage, the outlook for Planck is much better. Even without the highest-frequency chan- 
nels, which act as dust monitors, there is still a sufficiently broad spectrum of frequencies 
> 100 GHz where the spectral dependences of the DQ and the CMB quadrupole differ the 
most. Again, the synchrotron and bremsstrahlung foregrounds contribute mostly at low 
frequencies and thus do not degrade significantly the DQ signal. Quite remarkably, even if 
we marginalize over a number of uncertain foreground amplitudes, Planck should be able 
to detect the DQ as long as the velocity is greater than roughly 140 km sec -1 . Moreover, 
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TABLE I. Values for detectability of the dipole-induced intensity quadrupole. The affi column 
gives the standard error to the DQ amplitude. The \ 2 column quantifies the ability to distinguish 
v = from v = 627 km sec -1 using the DQ. The quadrupole components that are assumed 
to be fit to the multifrequency data are the DQ, CMB quadrupole (T), dust (D), synchrotron 
radiation (Synch), bremsstrahlung (Brem), far infrared background (FIRB), and Sunyaev-Zeldovich 
effect (SZ). The "NA" in the \ 2 column indicates that if the SZ effect is taken into account, our 
covariance matrix formally gives an infinite result because of the degeneracy between the frequency 
dependences of the CMB, the DQ, and Sunyaev-Zeldovich effect. The column Av gives the error 
(at 3a) to the magnitude of the velocity. We list NA for Av for the cases where the DQ is unlikely 
to be detected. The last column gives the anticipated (3a) error to the orientation of the velocity. 



experiment 


components 


a£Q (fiK) 


x 2 


Av (km sec x ) 


A9 (degrees) 


MAP 


DQ 


0.04 


10 5 


9 


0.9 


MAP 


DQ+T 


0.6 


440 


140 


14 


MAP 


DQ+T+Synch+Brem 


1.7 


52 


390 


42 


MAP 


DQ+T+D 


16 


0.66 


NA 


NA 


MAP 


DQ+T+D+Synch+Brem 


73 


0.032 


NA 


NA 


Planck 


DQ 


0.005 


6 x 10 6 


1.1 


0.12 


Planck 


DQ+T 


0.02 


3 x 10 5 


4.9 


0.52 


Planck 


DQ+T+D 


0.02 


3 x 10 5 


5.0 


0.53 


Planck 


DQ+T+D+FIRB 


0.036 


10 5 


8.1 


0.86 


Planck 


DQ+T+D+FIRB+SZ 


NA 


A 


NA 


NA 


Planck 


DQ+T+D+FIRB+Synch+Brem 


0.2 


4,000 


45 


4.8 
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Planck should be able to detect a deviation from the dip ole- inferred velocity as small as 45 
km sec -1 and determine its direction (modulo the sign) to better than 5°. 

It turns out that the Sunyaev-Zeldovich frequency dependence is a linear combination 
of the CMB-quadrupole and DQ frequency dependences. Thus, the DQ is formally degen- 
erate with the cosmological and Sunyaev-Zeldovich quadrupoles. However, the quadrupole 
moment due to the Sunyaev-Zeldovich effect will probably be negligible. Calculations of 
the Sunyaev-Zeldovich power spectrum near I ~ 100 find 1 2 Ci/(2tt) < 10 -12 |T7[. Doing 
shot-noise extrapolation to the quadrupole, we find cx^ z < 0.15 /xK, corresponding to a ve- 
locity uncertainty (at 3a) < 35 km sec -1 , which is just a bit below the expected statistical 
uncertainty indicated in Table | for Planck. The actual SZ number will probably be much 
smaller, and certainly much more will be known by the time Planck flies. The SZ octupole 
will be useful in constraining the possible SZ quadrupole contribution. 

Our foreground modeling makes several simplifying assumptions. For example, we do 
not actually know the frequency dependences perfectly and they can vary spatially (due to, 
e.g., spatial variation in chemical composition of the dust). There may also be components 
we have not yet considered, such as spinning dust. Further, about 25% of the sky will be 
lost to a Galactic-plane cut. It is thus possible that Planck may not be able to achieve 
the tiny velocity errors quoted in the table. However, the values in the final row may be 
achievable. Although uncertain, the frequency dependences of the foregrounds are all (with 
the exception of the SZ effect) considerably different than those of the the DQ and the CMB 
quadrupole. The foregrounds are all likely to have low amplitudes near 100 GHz: ~ 3 fiK 
for dust, ~ 1 fiK for bremsstrahlung and less for everything else ||18|| . Thus the DQ is not 
far below foreground quadrupoles, so multi-frequency foreground subtraction need not be 
done to better than about 10%. Of course, only the measurements themselves will answer 
these questions definitively. 



IV. AN INTRINSIC DIPOLE? 

Although a peculiar velocity is the simplest explanation for the dipole, and certainly that 
most consistent with the prevailing inflationary paradigm, it is also possible that the dipole 
could be due, at least in part, to an intrinsic temperature dipole produced by a super-horizon 
entropy perturbation ||. In this Section, we show that even if the temperature dipole is 
intrinsic, the same DQ still arises. We also clarify the distinction between an intrinsic CMB 
temperature fluctuation, the thermal quadrupole induced by a peculiar velocity, and the 
DQ, which unlike the other two, is not a temperature quadrupole. 

A velocity v induces a temperature pattern T{6) = Tcmb(1 + /3/i) -1 -\/l ~~ P 2 [0- To 
second order in /3 2 , this can be written as 

T(6) = T CMB {[1 + 0{(3 2 )} + /?V " 1/3)}- (5) 

Thus, the peculiar velocity induces a temperature quadrupole of magnitude 0(/? 2 T CMB ), in 
addition to the dipole of amplitude 0(/3T CMB ). However, the quadrupole in the intensity 
[Eq. (|3|)] arises from a combination of both the temperature dipole and the temperature 
quadrupole. To see this, we rewrite Eq. (||) as 
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T{6) = T CMB {[1 + 0{(3 2 )} - + - 1/3)}, (6) 
so that we can see where the DQ comes from. Doing so, we find that Eq. (|3]) becomes 

3 

^ = { I 1 + °^ 2 )] - f^f 3 ^ + + # M x ) - 1)] (/i 2 - 1/3) + •••}. (7) 

We thus see that what we have been calling the DQ consists of two parts: The first is the ft 2 
term, which is due to the temperature quadrupole induced by our peculiar velocity. Since this 
is an honest-to-goodness temperature quadrupole, it has the frequency dependence of the 
usual lowest-order thermal fluctuation, and it cannot be distinguished from an intrinsic CMB 
temperature fluctuation. The second term, proportional to (3f, arises from the term in the 
Taylor expansion of the intensity that is second order in the dipole amplitude. This second- 
order term has a frequency dependence that differs from the usual lowest-order thermal 
fluctuation. 

It is also clear from Eqs. (^) and (|7]) that even if the temperature dipole were intrinsic 
(that is, due to an entropy perturbation), then there would still be a DQ. Although the 
frequency dependence is g(x) — 1 rather than g(x), only the part proportional to g(x) can 
be distinguished by multi-frequency maps from a temperature quadrupole. For this reason, 
the amplitude and orientation of the DQ would be exactly the same as if the dipole were 
due to a velocity, and the detectability would be exactly as determined above. Likewise, 
the DQ cannot be used to tell whether the dipole is due to a peculiar velocity or due to an 
entropy perturbation, as previously suggested |12|.B 

Finally, we mention one last point. Strictly speaking, an intrinsic CMB temperature 
quadrupole will produce a contribution to the intensity quadrupole with a frequency de- 
pendence proportional to g(x), just like the DQ. However, if the temperature-quadrupole 
amplitude is AT (~ 10~ 5 for the intrinsic fluctuation and ~ 10~ 6 for the velocity-induced 
temperature quadrupole), then this contribution will be of order (AT) 2 < 10~ 10 , much 
smaller than the 10~ 6 expected for the DQ. 



V. CONCLUSIONS 

We have shown that the dipole amplitude is sufficiently large that the discrepancy be- 
tween the exact frequency dependence of a thermal fluctuation and the lowest-order fre- 
quency dependence usually assumed will be detectable by Planck and possibly MAP. To 
second-order in the dipole amplitude, this discrepancy is manifest as an intensity quadrupole 
that can be distinguished with multi-frequency measurements from an intrinsic temperature 
quadrupole. This provides a robust pre-determined target for CMB experiments and it may 
prove to be a useful tool for calibration for forthcoming space experiments. 

We thank the Santa Barbara KITP for hospitality. This work was supported at Caltech 
by NSF AST-0096023, NASA NAG5-9821, and DoE DE-FG03-92-ER40701, at Davis by 
NASA NAG5-11098, and at the KITP by NSF PHY99-07949. 
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